On the extinction of Continuous State Branching Processes 

with catastrophes 



Vincent Bansaye * 1 , Juan Carlos Pardo Millan t 2 , and Charline Smadi * 1,3 

1 CMAP, Ecole Polytechnique, Route de Saclay, F-91128 Palaiseau Cedex, France 
2 CIMAT A.C. Calle Jalisco s/n. CP. 36240, Guanajuato, Mexico 
3 CERMICS, Universite Paris-Est, 6-8 av. Blaise Pascal, Champs-sur-Marne, 

77455 Marne La Vallee, France 

January 28, 2013 



Abstract 

We consider continuous state branching processes (CSBP's for short) with ad- 
ditional multiplicative jumps, which we call catastrophes. Informally speaking, the 
dynamics of the CSBP is perturbed by independent random catastrophes which cause 
negative (or positive) jumps to the original process. These jumps are described by 
a Levy process with paths of bounded variation. Conditionally on these jumps, the 
process still enjoys the branching property. 

We construct this class of processes as the unique solution of a SDE and characterize 
their Laplace exponent as the solution of a backward ODE. We can then study their 
asymptotic behavior and establish whether the process becomes extinct. For a class of 
processes for which extinction and absorption coincide (including the a-stable CSBP 
plus a drift), we determine the speed of extinction of the process. Then, three subcrit- 
ical regimes appear, as in the case for branching processes in random environments. 
To prove this, we study the asymptotic behavior of a certain divergent exponential 
functional of Levy processes. Finally, we apply these results to a cell infection model, 
which was a motivation for considering such CSBP's with catastrophes. 
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1 Introduction 

Continuous state branching processes (or CSBP's for simplicity) are the analogues of 
Bienayme-Galton- Watson processes in continuous time and continuous state space. Such 
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classes of processes have been introduced by Jirina |Jir58j and studied by many authors 
included Bingham |Bin76] . Grey [Gre74| . Grimvall |Gri74j . Lamperti |Lam 67a, La m67bj . 
to name but a few. A continuous state branching process Z = (Zt,t > 0) is a strong 
Markov process taking values in [0, oo], where and oo are two absorbing states and 
(P x , x > 0) is the law of the process starting from x. Moreover, Z satisfies the branching 
property; that is to say, for any x, y > 0, P x + y is equal in law to the convolution of P x 
and F y . Since the pioneering work of Lamperti |Lam67b] . it is known that continuous- 
state branching processes are the only possible scaling limits of Bienayme-Galton- Watson 
branching processes and that every CSBP can be realized in this way. Thus CSBP may 
be model for the evolution of (renormalized) large populations which evolve during a 
large time window (see for instance [BTllj for parasite infection). 

The branching property implies that the Laplace transform of Zt is of the form 



Ex 



exp(— XZt) = exp{— xut(X)}, for A > 0, (1) 



for some non negative function ut- According to Silverstein [Sil68j . this function is deter- 
mined by the integral equation 

A 1 

— -du = t, (2) 

u t (\) VW 

where ip satisfies the celebrated Levy-Khintchine formula 

ip(X) = a\ + a 2 X 2 + / (e~ Xx - 1 + Xxl {x<1} ) fj,(dx), 

where a £ R, a > and u is a cr-finite measure on (0,oo) such that J( oo) \S ^ x 2 }n{dx) 
is finite. The function tp is known as the branching mechanism of Z. 

It is important to note that the first moment of Zt can be obtained by differentiating 
([T]) with respect to A. More precisely, note that d\Ut{X)\\ = Q = e - ^'( 0+ )* and E x (Z t ) = 
xe~^ '( 0+ )'. Hence, in respective order, a CSBP is called supercritical, critical or subcritical 
depending on ^'(0+) < 0, V'(0 + ) = or ^'(0+) > 0. Adding that 

lim Z t = 0] = e~ vx , 

where r\ is the largest root of the branching mechanism tp, the sign of ip'(0+) yields the 
criterion for a.s. extinction. Finally, we know that a CSBP Z with branching mechanism 
tp is absorbed at in a finite time with positive probability if and only if 

00 du 

< oo. 



In this case, P x (3t > : Z% = 0) = ¥ x (linit-j-oo Zt = 0) = exp(— rjx). 

In this paper, we focus on the case when Z has finite expectation, which is equivalent 

to 

/oo 
x/j,(dx) — a < oo, 

or equivalent to f^°(x A x 2 )n(dx) < oo. Then we can write 

^(A) = -gX + a 2 X 2 + J (e~ Xx - 1 + Ax) fi{dx), 
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and the CSBP is characterized by the triplet (g,a,fj>). Moreover a CSBP can also be 
defined as the unique non-negative strong solution of a SDE, which will be useful here. 
More precisely, from |FL10| . 

pt fit ft rco rZ a ^ ^ 

Z t = Z + gZ s ds + V ^Z s dB s + / / / zN {ds,dz,du), (3) 
Jo Jo Jo Jo Jo 

where B = (Bt, t > 0) is a standard Brownian motion, iVo(ds, dz, du) is a Poisson random 
measure with intensity dsfj,(dz)du independent of B and iVo is the compensated measure 
of N . 

The stable case is of special interest and will be studied in more detail in Section 3. It 
is motivated by applications (see Section 4) and provide us the key step to more general 
cases. It corresponds to the (only possible) scaling limit of Galton Watson processes 
(GW) when the reproduction law is fixed. Then 4>(X) = g\ + c\ a and two natural classes 
appear: 

• If a = 2 (Feller diffusion), we necessarily have fi = 0. Thus we get the continuous 
diffusion t 

Z t = Z + f gZ s ds + [ V^ 2 Z s dB s . 
Jo Jo 

• If a S (1)2), we have that a = and /j,(dx) = cax~ a+1 dx/T(2 — a). Then the 
process has (only) positive jumps with infinite intensity |Lam07| . Moreover, we 
have ^ ^ 

Z t = Z + I gZ s ds+ [ z]L a dX s , 
Jo Jo 

where X is an a-stable spectrally positive Levy process. 

In the stable case, we can solve equation ([2]). In particular, we know precisely the speed 
of absorption, i.e. the asymptotic behavior of F(Z t > 0). 

In this work, we want to take into account catastrophes which occur randomly and 
kill each individual with some probability (depending on the catastrophe). In the scaling 
limit of large population (the continuous state setting), it amounts to let the process 
make a negative jump and multiply Z by a random fraction. The process we obtain is 
still Markovian if the catastrophes happen without memory and the fractions are i.i.d., 
in other words catastrophes are described by a Poisson Point Process. Conditionally on 
the times and the effects of the catastrophes, the process Z also satisfies the branching 
property, as detailed below. It yields then a particular class of CSBP's in random envi- 
ronment, which can also be obtained as scaling limit of GW in random environment (see 

[EM]). 

More precisely, we consider here CSBP's with catastrophes described by an independent 
Poisson Point Process. Such processes are motivated in particular by cell division mod- 
els. For example, in [BTllj . the case of Feller diffusion with catastrophes describes the 
evolution of parasites in a cell line. The parasites grow following a Feller diffusion and 
the division of the cell causes a negative jump of this quantity by splitting the parasites 
into two parts. 

More generally, we can also take into account positive jumps due to favorable environ- 
ments. They can represent immigration events proportional to the size of the population. 
One possible application comes from the aggregation behavior of some species. Indeed, 
many species form aggregates, what may allow them to reduce predation, make the food 
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search easier, or improve the mate choice (see respectively chapters 12.1, 12.2.2, 12.2.3 
and 12.3.4 in |DGU08| ). For convenience, we still call the processes obtained CSBP with 
catastrophes. 

The remainder of the paper is structured as follows. In the next section, we study the 
CSBP with catastrophes, in particular we verify that these kinds of processes are well 
defined by a SDE. We also characterize their Laplace exponent via an ODE which allows 
us to describe their long time behavior. In particular, we determine when (Yt,t > 0), a 
CSBP with catastrophes becomes extinct. That is, we specify when P(limt_ s . 0O Yi = 0) = 1 
and the asymptotic behavior when the process survives, under some moment assumption. 

Section [3] is devoted to the speed of extinction. We take special attention to the 
stable case (Section 13. ip . where the extinction event coincides with the absorption event. 
In particular, we obtain an explicit expression for the Laplace exponent in terms of an 
exponential functional of a Levy process whose jump structure is given by the catastro- 
phes. This expression allows us to determine the speed of absorption of the process when 
it is subcritical or critical (see Section [5]). 

From this result in the stable case, we can deduce the speed of extinction for a large 
class of CSBP's with catastrophes, when extinction and absorption coincide. In the 
subcritical case, three regimes appear (see Theorem [S] and Corollary . It means that 
the law of the catastrophes (and not only their mean effect) can change the exponential 
decrease of the survival probability. This result is closely related to the discrete framework 
via asymptotic behavior of functionals of random walks. More precisely, we use in our 
arguments local limit theorems for semi direct product |LPP97llGL01j and some analytical 
results on random walks |Koz76t IHir98j (see Section [5|). In the same vein, we refer to 
|BH12j for different subcritical regimes for Feller diffusions which evolve "in Brownian 
environment" . 

This result can be related to discussions in ecology about the role of environmental 
and demographical stochasticity. Such topics are fundamental in ecology and conservation 
biology, as discussed for instance in Chapter 1 in |LES03j . In our model, the survival may 
be either due to the randomness of the individual reproduction, which is given by the 
parameters a and fi of the CSBP, or to the randomness (rate, size) of the catastrophes 
which is linked to the environment. For a study of relative effects of environmental and 
demographical stochasticity, one can read |Lan93| and references therein. 

In Section 01 we apply our results to Feller diffusion with catastrophes and derive the 
different regimes for the speed of propagation of an infection in the cell division model 



Section 5 is devoted to the asymptotic behavior of exponential functionals of Levy 
processes with paths of bounded variation. More precisely, we are interested in the 
asymptotic behavior at oo of 



where /3 > 0, F belongs to a given class of functions (see (|26p ) and K is a Levy process 
with paths of bounded variation that does not drift to oo and has some exponential 
positive moments. We find four different regimes for its behavior that depend on the 
shape of the Laplace exponent of the process K. Up to our knowledge, this situation has 
been studied only by Carmona et al. |CPY97j . Lemma 4.7, where they provide only one 
precise regime. The case when the process K drifts to oo has been deeply studied by 
many authors, see for instance Bertoin and Yor |BY05| and references therein. 



[BTllj . 
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Finally, Section 6 contains some technical results which are used in the proofs. We 
establish such results at the end of the paper for the convenience of the reader. 



2 CSBP with catastrophes 

We consider a CSBP Z = (Z t ,t > 0) defined by Q characterized by the triplet 
(g, a, fi), where we recall that [i satisfies 



oo 



(x A x 2 )n(dx) < oo. (4) 







The catastrophes are independent of the process Z and are given by a Poisson Point 
Process (5 eti ,ti)ie/ on [0, oo) x [0,oo) with intensity dtv(dx) such that 

(1 A | a; - l|)z/(d£c) < oo. (5) 

(0,oo) 

The associated Poisson random measure Ni = X^iei^et-,*, is independent of B and No, 
and the jump process 

A t = f f log(x)iVi(ds,dx) = Vlog(e s ), 

JO J(0,oo) s < t 

is a Levy process with paths of bounded variation. 

The CSBP (g, a, fi) with catastrophes v is then given by the following SDE 



Y t = Y + [ gY s ds+ [ V^ 2 Y s dB s + f [ [* zN {ds, dz, du) 

Jo Jo Jo J[0,oo) Jo 

(z-ljYs-N^dz), (6) 



J[0,oo) 

where Yq > a.s. In other words, when a catastrophe occurs at time s, the population 
size is multiplied by e s which is distributed according to the measure v. 

Theorem 1. The SDE has a unique non-negative strong solution for g £ R, a > 0, 
fj, and v satisfying the conditions ^> and respectively. 

Then, the process Y = (Yt,t > 0) defined by the SDE @) is a cadlag Markov process 
satisfying the branching property conditionally on A = (At,i > 0) and with generator 
given by 

Af(x) = gxf'(x) + <J 2 xf"(x) + (f(xz) - f(x)) u{dz) 

J ° V 7 N (7) 

f(x + z)-f(x)-zf'(x))xn(dz). 



+ 

Moreover, for every t > 0, 



o 



Ey 



expj - Aexp{ - gt - A t }lt| A = exp j - yv t (0, A, A)| a.s 
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where for every (A, 5), v t : s t— > vt(s, A, S) is the unique solution of the following backward 
differential equation : 

^- s V t (s,X,5)=e3 s+5 ^ (e-3 s -^v(s,\,5)), 0<s<t, v t (t,X,5) = X, (8) 

where 

roc 

V>o(A) = ^(A) - A^'(O) = a 2 A 2 + / (e~ Xx -l + \x)n(dx). (9) 



Jo 

Proof. We first study the SDE given by 

-Y. 



Y t = Y + [ gY s ds+ [ V^ 2 Y s dB s + [ [ f' zN (ds,dz,du) 

JO JO JO J[0,oa) Jo 

(z-ljYs-N^ds^z). (10) 



J (0,2) 



The above SDE satisfies the conditions of Theorems 3.2 and 5.1 of |FL10j . In particular, 
to verify condition (5.b) we can simply let V n = [1/ra, oo) x (0, oo), for n > 1. Hence 
from Theorems 3.2 and 5.1 in [FLIOj . we have pathwise uniqueness and the existence of 
a unique non-negative strong solution for (|10jh 

Now, note from © that v([2, oo)) < oo. Thus, Proposition 2.2 in [FLIP] implies pathwise 
uniqueness and the existence of a unique non-negative strong solution for the SDE ([6]). 
By ito's formula (see for instance |IW89j Th.5.1) one sees that the solution (Yt,t > 0) 
solves the following martingale problem: for every / £ C 2 (IR+), 

f(Y t )=f(Y )+ loc. mart. + g f f'(Y s )Y s ds 

Jo 

pi /""t p oo 

+ a 2 J f"{Y s )Y s ds + J J Y s [f{Y s + z) - f(Y s ) - /'(y s )z)/i(dz)d S 
f{zY s -)- f{Y s -))u(dz)ds, 



t fOO 



+ 

'0 JO 

where the local martingale is given by 

ct ft POO 



PL PL OO 

loc. mart. = J f {Y s )^2oW s dB s + J J (f(zY,-)-f(y a -j)N 1 (ds,dz) (11) 

pi /* oo /* y 

+ J J J S (f(Y s _ + z)-f(Y s ^N (ds,dz,du), 



and N\ is the compensated measure of N%. From pathwise uniqueness, we deduce that 
the solution of © is a strong Markov process whose generator is given by ([7J). 

The branching property of Y, conditionally on the jumps, is inherited from the branch- 
ing property of the CSBP and the fact that the additionnal jumps are multiplicative. 

In order to prove the second part of the theorem, let us now work conditionally on A. 
By applying Ito's formula to the process Zt = Yte~ gt ~ At , Lemma [TBI in Section [6] ensures 
that for every F £ C^' 2 , F(t, Zt) is also a local martingale if and only if for every t > 0, 

pt £l2 pt £J 

j^—F(s,Z s )a 2 Z s e^ s - A °ds + J -F(s,Z s )ds 

pi P OO Q 

+ J J Z s (^F(s,Z s + ze-° s - A °) - F(s,Z s )]eV s+A ° - — F(s,Z s )z^(dz)ds = 0. 
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In the vein of |IW89[ [BTllj . we choose F(s,x) := exp{— xv t (s, A, A)}, where vt(s,X,A) 
is differentiable with respect to the variable s, non negative and such that vt(t, A, A) = A, 
for A > 0. We observe that F is bounded. Therefore, from the above computations 
(exp{— Z s vt(s, A, A)}, < s < t) will be a martingale if and only if 



d , \ \ \ 
—v t (s,X,A) = 



gs+A t 



■00 (e 



-gs-A. 



v t (s,X,A)), 0<s<t, 



where ipo is defined in ([9]). Proposition [16] in Section 6 ensures that a.s. the solution of this 
backward ODE exists and is unique, which essentially comes from local Lipschitz property 
of ipo and paths of bounded variation of A. Thus the process (exp{— Z s vt(s, A, A)}, < 
s < t) is a martingale conditionally on A. We get 



E, 



which implies 



exp | - Z t v t (t, A, A) | 



A 



E, 



exp | - Z o v t (0,X,A)j 



A 



E, 



exp|-AZ t | A = exp | -yv t (0,X, A)|, 



and ends up the proof. 



(12) 
□ 



Now we derive from Theorem [T] an extinction criterion for CSBP's with catastro- 
phes. First, we recall that a Levy process L = (Lt,t > 0) has only three types of 
asymptotic behaviors: either it drifts to oo, — oo, or oscillates a.s. The latter means that 
limsup^oo L% = — lim inf t-Kxj Lt = oo. We refer to Theorem 7.2 in |Kyp06| for a criterion 
for such types of behaviors. 



Corollary 2. We have the three following regimes. 



i) If {A t + gt)t>o drifts to — oo, then Yf — > a.s. 

ii) If (At + gt)t>o oscillates, then liminft_s.oo Yt = a.s. 

Hi) If (At + gt)t>o drifts to +oo and there exists e > 0, such that 

POD 

/ xlog(l + x) 1+e n(dx) < oo, (13) 
Jo 

then P(liminf i _ i . 00 Yt > 0) > and there exists a non negative finite r.v. W such that 
e -gt-A tYt >w a ^ w = Q | = | Hm y t = o). 

It is important to note, according to Theorem 7.2 in |Kyp06| , that whenever Ai has a 
finite first moment, then the sign of E(Ai) +g yields the regimes above. Hence, whenever 
its Laplace exponent, (j)(X) = logE[exp(AAi)], is well defined for some positive A, then 
cases (i), (ii) and (Hi) will be equivalent to g + ^'(0) < 0, g + (j)'(0) = and g + ^'(0) > 
respectively. This will be the case in the rest of the paper. 

Moreover, let us note that in the regime (ii), we can have lim inf t->oo Yt = a.s. but 
Yt a.s. does not tend to zero. For example, if \i = and a = 0, Yt = exp(gt + At), hence 
limsup^oo Yt = oo. 



7 



Assumption (Hi) of the corollary does not imply that {lim^oo Yt = 0} = {3t : Yt = 0}. 
For example if /i(dx) = x~ 2 l[ ^dx, it is not difficult to show that ip(u) ~ ulogu as 
u —> oo, and according to Remark 2.2 in |Lam08| we necessarily have P(3i : Yt = 0) = 
< Pflimt-Hx, r t = 0) < 1. 

Proof. We know from ()12|) with -F(s, x) = x, that exp{— gt — At}lj is a non negative 
local martingale, thus it is a non negative supermartingale and it converges a.s. to a non 
negative finite random variable W. This leads to the result for the first two cases. 
In the case when (At + gt)t>o drifts to +oo, we should now prove that F(W > | A) > 
a.s. According to Lemma [181 in Section 6, there exists a non negative increasing function 
h on M + such that for all A > 0, 



V>o(A) < A/i(A) and c(A) := J /i(V (9 * +At) )d£ < 



oo a.s. 



We can now show that conditionnally on A there exists a positive lower bound for 
vt(0, 1, A), for t > 0. For every (t, A) G the solutions of ([8]) are non decreasing on 

[0,t], thus for all s G [0,t], u f (s, 1, A) < 1, and 

Me- 9S ~ As vt(s, 1, A)) < e-° s - As v t (s, 1, A)/ l ( e -^" As ^( S , 1, A)) 

< e^ As v t (s,l,^)He~ 9S ~ As ). 

Then 

^^(s, 1, A) < u t (s, 1, A)/i(e" 9S_As ). 

This clearly implies, 

- In(«t(0, 1, A)) < / /i(e- 9S - As )(is < c(A) < oo, a.s. 

JO 

Hence, for every t > 0, u t (0, 1, A) > e _c(A) > 0. Then from ([T2~]) . 

E 3/ (exp(-AVF) | A) = exp ( - y lim u t (0, 1, A)) < 1 



and P(PF > | A) > 0. 

Moreover, since Y satisfies the branching property conditionally on A, we can show 
(see Lemma [T9l in Section 6) that 

{W = 0} = f lim Y t = o\ a.s., 

which completes the proof. □ 

We now derive a central limit theorem in the supercritical regime that only requires 
a second moment assumption. 

Corollary 3. Assume that (At + gt)t>o drifts to +oo and i!3\) is satisfied. Then, under 
the additional assumption 

(log(a;)) 2 ^(dx) < oo, (14) 

(0,e- 1 ]u[e,oo) 



S 



conditionally on {W > 0}, 

log(y t ) - mt d 



pVi 



t—too 



> N(0, 1) 



where A means convergence in distribution, 

m:=g+ I log(x)i/(dx) < oo, p 2 := / (log(x)) 2 ^(dx) < oo, 

J{\log(x)\>l} Jo 

and N(0, 1) denotes a centered gaussian random variable with variance equals 1. 

Proof. We first establish a central limit theorem for the Levy process (gt + At,t > 0) 
under the assumption (|14p which relies on Theorem 3.5 in Doney and Mailer [DM02J (the 
details can be found in Section [6~5]) . We get 



^ + A '" mt d >iV(0,l). (15) 



From Corollary [2] part Hi), on the event {W > 0}, we then have 
log(r t ) - (gt + A t ) log(^) E (-00,00) 

t— >oo 

and we conclude by using (fT5|) . 



□ 



3 Speed of extinction of CSBP with catastrophes 

In this section, we first consider the stable CSBP with a deterministic growth g S R, 
before deriving the result for some general class of CSBP. 

3.1 The stable case 

We assume in this section that 

V(A) = < 7 A + c + A /3+1 , (16) 

for some f3 £ (0, 1], c_|_ > and g in BL In this particular case, the backward differential 
equation (JSJ) can be solved and we get 



t 

/ e -p(A s+gs ) ds 



Proposition 4. For all xq > and t > 0: 

F X0 (Y t >0) = 1-E 
Moreover, 

F Xo (there exists t > 0; Y t = 0) = 1, 
if and only if the process (At + gt,t > 0) (ioes not drift to +00. 



(17) 
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Proof. We solve equation (|SJ) with tp(\) = gX + c + A^ +1 . As Y>o(A) = c+A^ +1 , a direct 
integration yields 



v t (u,t, A) 



c+ f e-^ A * + ^ds + A" 

■/ u 



■p 



1-1/13 



so that 



E 



exp < -x 



c+ / e^^ds + A"' 3 



(18) 



The expression for the absorption probability is a direct application of ([T8]) . Indeed letting 
A goes to oo, we get 

e -/3(A s+5 s) ds 







PzoOt = 0) = E x . 
The mapping 1 1— >• J e~^^ As+9S ^ds is a.s. increasing. Thus, 

e -/3(A s +gs) ds G M +U { 00 } ) 



f 

JO 



is well defined. Hence, the process (Zt,t > 0) converges in distribution as t — >• oo towards 
the r.v. W (already defined in the case Hi) of the Corollary whose distribution is 
specified by 



E 



■CO 



-AW 



E 



exp I -xq ( c+ / e -ft*" +gs) ds + A' 



Letting A — > oo, we get by monotone convergence 
P a;o (W = 0) =E X0 



exp < -x c+ 



oo x -1//9 

e -/3(A s+5 s) ds 



(19) 



Finally, according to Theorem 1 in |BY05| we have J °° e ^( As +# s )d,s = oo, a.s. if and 
only if the process (At + gt,t > 0) does not drift to +oo. Adding that ¥ X0 (3t > 0; Yt = 
0) = \mx¥ xa (Y t = 0) and using flT7J) and flUD, we get that 

t— >00 

1 = F X0 (W = 0) = F, (3i > 0; Y t = 0) 

if and only if the process (At + gt,t > 0) does not drift to +oo. □ 

In what follows, we assume that the Laplace exponent of the Levy process A is well- 
defined for some positive real number, i.e. 

0(A) = log E[e AAl ] for A G [0, 9 max ), 

where 9 max = sup{A > 0, 0(A) < oo}. In other words, Jj e x x x v(dx) < oo if A G [0, 9 max ). 
We get then the main result of this paper, where we recall that (f)' is non decreasing. 



Theorem 5. We assume that ip satisfies Mb}) . 
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a/ If 4>'(0) + g < (subcritical case) and O max > 1, then we have the following three 
asymptotic regimes 

(i) If 4>'(1) + g < (strongly subcritical regime), then there exists c\ := 
di(/3, zv)^ 1 ^ > such that for every xq > 0, 

P X0 (Y t > 0) ~ Cl x e mi)+9 \ as t oo. 

(ii) If (f)'(l) + g = (intermediate subcritical regime), then there exists C2 := 
g?2(/3, ^)c + 1 ^ > such that for every xq > 0, 

F X0 (Y t > 0) ~ c 2 x t- 1 / 2 e t ^ +9 \ as t oo. 

(m,) If (p' (1) + g > (weakly subcritical regime) and 6 max > (3 + 1, i/ien /or every 
.xo > 0, t/iere exists C3 := €3(^0, V 7 )^) > suc/i t/iai 

F X0 (Y t > 0) ~ cgi" 3 / 2 ^^^, as f — »■ 00, 

where r is f/ie rooi of <p' + g on ]0, 1[. 

0/ If <j)'{0) + 5 = C0 (critical case) and 6 max > (3, then for every xq > 0, i/iere exists 
C4 := C4(xo,V') zy ) > suc/i t/iat 

^o^t > 0) ~ c 4 t" 1/2 , as t^co. 



Proof. From the Proposition [5] we know that 



p xo (y t > 0) = i-E 



exp 



-x 









l 


= E 











where F(x) = 1 — exp{— xq(c + x)~ 1 /^} and K s = A s + (7s. 

Since F satisfies (|26|) . the result is a consequence of Proposition [8] whose statement and 
proof are deferred to Section [5j □ 

In the case of the CSBP without any catastrophes, the subcritical regime is reduced 
to (i). The critical case differs from b/, since then the asymptotic behavior is given by 
1/t. 

This result can be compared to the regimes which appear in the literature of 
discrete (time and space) branching processes in random environment (see e.g. 
[GLOU IGKV03| IAGKV05] ). even if the proofs do not use directly the corresponding 
results. 

Besides, in the continuous framework, such asymptotic regimes for the survival prob- 
ability of Feller diffusion whose drift varies following a Brownian motion appear in |BH12j . 



In the strongly and intermediate subcritical cases (i) and (ii), E(y) gives the good 
exponential rate for the decrease of the survival probability, + g. Moreover the 

We exclude the degenerated case v = 0, g = 0. 
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probability of non-extinction is proportional to the initial state xo of the population, and 
to c, 1 ^ . This observation follows from the following equivalence 



E 



FU exp(-/3(A s+5 s)) 



E 



x c + 



[ exp(-/3(A s + 5 s)) 
J o 



as 



oo, 



and we refer to the proof of Lemma [IT] and Section [5.41 for details. 

In the weakly subcritical case (Hi), the exponential rate of decrease of the survival 
probability is <p(r) + gr. The latter is strictly smaller than 0(1) + g, indeed 

4>{t) + gr = mm {(/)(s) + gs} < 0(1) + g. 

0<s<l 

In fact, as it appears in the proof of this theorem (see Section [5]), the quantity which 
determines the scale of the asymptotic behavior in all cases is linked to E[e /j ], where 
It := mf sem (A s + gs). 

Let us note also that C3 and C4 may not depend linearly on xq. We refer to |Ban09j for a 
result in this vein for discrete branching processes in random environment. 



3.2 Beyond the stable case. 

In this section, we want a similar result as in Theorem [5] for CSBP's with catastrophes 
when the branching mechanism ipQ is not stable. In order to do so, we may compare ipo 
with a stable branching mechanism. For some technical reasons, we assume that the 
Brownian coefficient is positive and the associated Levy measure /i satisfies a second 
moment condition. It allows us to obtain the following result from the Feller diffusion 
case in Theorem [5l i.e j3 = 1. 

Corollary 6. Assume that J/ Q s x 2 [i(dx) < 00 and a 2 > 0. 

a/ In the subcritical case, i.e. 4>'(0) + g < 0, if 9 max > 1, we have the following three 
asymptotic regimes 

(i) If <p'(l) + g < 0, there exist < c\ < c[ < 00 such that for every x$, 

cix e* w(1)+9) < ¥ X0 (Y t > 0) < c[x e tWl)+9) (t > 0). 

(ii) If <p'(l) + g = 0, there exist < C2 < c' 2 < 00 such that for every xq, 

C2 x i- 1/2 e^«+^ < F X0 (Y t > 0) < c^ot^e^ 1 ^ (i > 0). 

(Hi) If 4>'(l) + g > and 6 ma x > P + 1> for every xq, there exist < 03(^0) < 
c' 3 (xq) < 00 such that 

c 3 (*o)i~ 3/ V {< « T)+5T) < F X0 (Y t > 0) < c' 3 (x )t- 3 / 2 e t ^ +9T '> {t > 0), 

where r is the root of <p' + g on ]0, 1[. 

b/ In the critical case, i.e. <j)'(0) + g = dl, if 6 max > ft, then for every x$, there exist 
< C4(xo) < 04(^0)' < co such that 

c 4 (x )t- 1/2 < F X0 (Y t > 0) < 4(x )t" 1/2 (t > 0). 
2 We exclude the degenerated case v = 0, g = 0. 
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We observe that in the proof, the assumption a 2 > is required only for the upper 
bounds in the previous inequalities. 

Proof. We recall that the branching mechanism associated with the CSBP Z satisfies, for 
every A > 

^(A) = -gX + a 2 X 2 + J (e~ Xx - 1 + Ax) fi(dx). 



So 

V>"(A) = 2a 2 + I x 2 e- Xx fi(dx), 

J(0,oo) 

and for every A > 

2a 2 < i/j"(X) < 2a 2 + f x 2 e- Xx fi(dx). 

J(0,oo) 



r (0,oo) 

Since c := J °° x 2 [i(dx) < oo, then tp" is continuous over [0,oo). By Taylor Lagrange 
theorem, we get for every A > 0, ip~(X) < ip(X) < ?/>+(A), where 

z/>_(A) = A^'(0) + a 2 X 2 and r/> + (A) = A^'(0) + (a 2 + c/2)X 2 . 

We first consider the case z/(0, oo) < oo, so that A has a finite number of jumps on 
each compact interval a.s. We introduce the CSBP's with catastrophes Y_ and Y + which 
are both associated with the same catastrophes as Y, A, but respectively with the CSBP 
(g, cx 2 ,0) and (g,a 2 + c/2,0). We denote U-j and their Laplace exponent, i.e. for all 
(A, t) G 

E exp{-AY t ~} = exp{-tt_ t (A)}, E exp{-AY t + } = exp{-u +it (A)}. 

Thus conditionally on A, for every time t such that At = A^_, by Theorem [TJ 

u '-,tW = -^-( u -,t)i u '+,tW = u' t (X) = -ij)(ut). 

Moreover for every t such that 6t = expjA^ — A^_} > 0, it is clear that 

u-A X ) = u tW = U +A X ) = e 
u- )t -{\) tit- (A) u+,t-(A) 

and U-fi{X) = uo(X) = -u +i o(A) = A. So for all t, A, 

u+,t(A) < u(t,X) < u_ lt (A). 

Now we generalize this inequality to the case z/(0, oo) G [0, oo] by successive approxi- 
mations. With this purpose, let A ei,€2 = (0, 1 — s\) U (1 + £2, oo), where 0<1 — ei < 1 < 
l+£2 and define the Poisson random measure N^ 1 ' 62 as the restriction of N\ to A £1 ' 62 x M + . 
We denote by dtu £l,£2 (dx) its intensity measure, where z^ £l ' £2 (dx) = l{ X £An,£2}v(dx), and 
the corresponding Levy process A £l,£2 defined by 



A euS2 = f* I" l g(x)iV 1 ei ' e2 (ds,dx). 
JO J(0,oo) 



'(0,oo) 

We also consider the CSBP's Y £l ' £2 (resp Y £l ' £2 ~ and Y £l ' £2 ' + ) with branching mech- 
anism ip (resp. tp- and and catastrophes A £l,£2 via Since z^ £l,£2 (0, oo) < oo, 
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from the first step we have u^:f 2 (X) < u £l ' 62 (t, X) < nl 1 ' i £2 (A), where as expected 
Efexpl-AY/ 1 ' 22 '*}] = exp{-<;f 2 '(A)} for each * G {+,0,-}. ' 

Similarly, let A ei = (0, 1 — ei) U (l,oo) and define the Poisson random measure N^ 1 
as the restriction of TVi to A ei x IR + whose intensity measure is given by dtv 6l (dx), where 
v 61 (dx) = li x eA e i} u {dx). Let us fix t in W* + , and define Y 61 the unique strong solution 
of 

Y t £1 =Y + [ gY s £l ds+ [ y/2(T 2 Y 8 ei dB a + I f [ " zN {ds, dz, du) 

Jo Jo Jo J[0,oo) Jo ^Q) 

+ 11 (z-\)Y^N^[ds,dz). 

Jo J[0,oc) V J 

We already know from Theorem [T] that equation (|2U|) has a unique non negative strong 
solution. Moreover, from Theorem 5.5 in [FLIP] and the fact that N 61 has the same 
jumps as N El ' e2 plus additional jumps greater than one, we know that 

y/ 1,£2 <Y/\ a.s. 

Using assumption (|5|), we can apply Gronwall Lemma to the non negative function t i— > 

E[Y t £l - Y t ei ' E2 \ and obtain 



E 



£2^-0 



■> 0. 



Adding that Y £l ' E2 is decreasing with £2, we finally get, Y^ 1,e2 Yf 1 , as £2 - > 0. 
Using similar arguments as above for Y ei ' e2 ' + and Y £l ' £2 ~, we deduce 

u%(X)<u^(t,X)<u%(X). 

Thus, when e\ goes to 0, we finally obtain u+ t t(X) < u(t,X) < u_^(A). This implies, 
taking A — > 00, that 

F(Y t + > 0) < F(Y t > 0) < F(Y f - > 0). 
The result now follows from the asymptotic behavior of F(Y t ~ > 0) and F(Y t + > 0). □ 



4 Application to a cell division model 

When the reproduction law has a finite second moment, the scaling limit of the Galton 
Watson process is a Feller diffusion with growth g and diffusion part a 2 . It yields the 
stable case with /3 = 1 and additional drift term g. Such process is also the scaling 
limit of birth and death process. It gives a natural model for populations which die and 
multiply fast, randomly, without interaction. Such a model is considered in [BTllj for 
parasites growing in (dividing) cells. In this model, the cell divides at constant rate r 
and a random fraction in (0, 1) of parasites goes in the first daughter cell, whereas the 
rest goes in the second daughter cell. Following the infection in a cell line, the parasites 
grow as a Feller diffusion process and undergo a catastrophe (at finite rate) when the cell 
divides. If there is one infected cell at time 0, the numbers Nt of cells and N£ of infected 
cells at time t satisfy E[N t ] = e rt and E[JV t *] = e rt F(Y t > 0), where 

Y t =l+ [ gY s ds+ [ ^2a^Y s dB s + [ [ (8 - l)Y s _p(ds,d6). (21) 
Jo Jo Jo Jo 

Here B is a Brownian motion and p(ds,d0) a Poisson Point measure with intensity 
2rdsF(Q G d6). Then Theorem [5] and Corollary [5] directly ensure the following result. 
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Corollary 7. a/ We assume that g < 2rE [log(l/0)] . Then there exist positive con- 
stants c\ , C2 , C3 such that 

(i) Ifg< 2rE[eiog(l/0)], then 

E [Nf] ~ c x e 9 \ as t -> 00. 
// 9 = 2rE[01og(l/e)], ffcen 

E [iV t *] ~ c 2 t" 1/2 e 9 *, as i -> 00. 
(mj Ifg > 2rE[01og(l/0)], then 

E [JV t *] ~ c 3 t" 3/2 e at , as i -> 00. 
w/iere a = min Ae [ ,i]{fl^ + 2r(E[0 A ] - 1/2)} < g. 
b/ We assume now g = 2rE [log(l/0)], then there exists C4 > such that, 

E [JV t *] ~ c 4 t" 1/2 e r *, as t -> 00. 

c/ Finally, if g > 2rE [log(l/0)], then there exists < C5 < 1 such that, 

E[iV t *] ~ c 5 e rt , as i -> 00. 

Hence if 5 > 2rE [log (1/0)] (supercritical case c/), the mean number of infected cells 
is equivalent to exp(ri), which is the mean number of cells. In the critical case (b/), there 
are a bit less infected cells, owing to the additional square root term. In the strongly 
subcritical regime (a/ (i)), the mean number of infected cells is of the same order as the 
number of parasites. It let us think that parasites do not accumulate in some infected 
cells. The asymptotic behavior in the two remaining cases is more complex. 

We stress on the fact that fixing the growth rate g of parasites but making the law of 
the repartition vary make change the asymptotic behavior of the number of infected 
cells (when E[01og(l/0)] goes beyond g/2r). For example, let us fix the growth rate 
of parasites g and the rate of division r of the cells such that g < rlog2. It yields 
g < 2rE [0 log(l/0)] for = 1/2 a.s. We make now increase the entropy of the repartition 
of parasites E [0 log(l/0)] (the function x 1— >■ x ln(l/x) is strictly concave on [0, 1] and null 
at the boundaries). In the domain g < 2rE [01og(l/0)], the growth rate of the number 
of infected cells is unchanged. A threshold appears when E[01og(l/0)] = g/2r. Then, 
this growth rate decreases and the number of infected cells grows slower and slower. Such 
phenomena have already been observed in the discrete time, discrete space framework in 
[Ban08| . 

The results can also be extend to a growth of parasites which follows a stable CSBP. 
Such model appears when renormalizing a discrete model where the reproduction law of 
parasites has an heavy tail. 

5 Local limit theorem for some functionals of Levy pro- 
cesses 

We consider a Levy process K = (Kt, t > 0) of the form 

K t = -ft + cjJ +) - <r\~ ] , t>0, (22) 
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where 7 is a real constant, a^ + > and > are two independent subordinators without 
drift. We denote by II, IlW and II^ - ) for the associated Levy measures of K, <r(+) and 



a 



(-) 



Our aim, in this section, is to determine the asymptotic behavior of the distribution 
of the exponential functional associated to K, 



/ exp{-/3K s }ds, 
Jo 



where j3 belongs to (0, 1]. 

Let us define the Laplace exponents of K, cj( + ) and cr(~ ) by 



<f> K (\) = logE[e AKl ], 0+(A) = logE 
and assume that 



(+) 



and 4>- (A) = log E 



fmax = SUp 

From the Levy-Khintchine formula, we deduce 



JagM + , / e Al n( + )(dx) < 00 I > 0. 
{ J [hoo) J 



<^(A)= 7 A+ / (e Xx -l)u^(dx)+ [ (e~ Xx -l)rt~\dx) 

J(0,oo) V ' J (0,oo) ^ ' 



'(0,oo) v ' J (0,oo) 

Finally, we assume that K(Kf) < 00, which is equivalent to 

x 2 Yi{dx) < 00. 



00,00) 



(23) 



(24) 



(25) 



We are interested in the asymptotic behavior at infinity of 



a F (t) :=E 



F (J exp{-PK„}ds 



where F belongs to a particular class of functions on R_|_. We will focus on functions which 
decrease polynomially at infinity (with exponent —1//3) in order to get a fine estimate 
which yields the four asymptotic regimes in Theorem [5J 

Proposition 8. Assume that F is a positive non increasing function such that for x > 



F(x) = C F {x + 1)~ 1//3 1 + (1 + x^hix) 



(26) 



where ? > 1, Cp is a positive constant, and h is a Lipschitz function which is bounded. 
Then we have the four following regimes 

a/ (Subcritical case) If 4>' K (0) < 

(i) If 9 max > 1 and 4>' K (l) < 0, then there exists a positive constant c\ such that 

a F (t) ~ cie ife(1) , as t -> 00. 
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(ii) If 6 max > 1 and </>/£-(!) = 0, then there exists a positive constant C2 such that 

a F (t) ~ c 2 t" 1/2 e*^' (1) , as t ^ oo. 

(Hi) If 6 max > + 1 a^d ^(1) > 0, t/ien t/iere exists a positive constant C3 suc/i 
that 

a F (t) ~ c 3 t- 3/2 e^ K{r) , as t -> 00, 
where r is i/ie rooi of 4>' K on ]0, 1[. 

6/ (Critical case) If Q max > P &nd </>^(0) = ; £/ien i/iere exists a positive constant C4 
suc/i i/wi 

a_F(i) ~ C4t -1 ^ 2 , as t — > 00. 

This result is proved in Section 15.41 It generalizes Lemma 4.7 in Carmona et al. 
[CPY97] in the case when the process K has paths of bounded variation. We also note 
that the result of Carmona et al. does not have a precise asymptotic behavior when 

The behavior of a^(t) is determined by the shape of (ftx, which is the Laplace exponent 
of K. This function is convex and zero at the origin and four different regimes appear 
according to the sign of 4>' K on [0, 1]. 

The assumption on the tail of F as x — > 00 is finely used to get the asymptotic 
behavior of ap (t) . Lemma [21] gives the properties of F which are required in the proof. 

The strongly subcritical case (case (i)) is proven using a continuous time change of 
measure. For the other cases, we divide the proof in three steps. The first step consists in 
discretizing the integral e~^ Ke ds by introducing for every (p, q) £ N x N* the random 
variable 

p p i—l 

A m = exp{-pK i/q } = ^[]exp{- P{K {]+l)/q - K j/q )}. (27) 

j=0 i=0 j=0 

Secondly, we study the asymptotic behavior of the discretized expectation : 

F(A p Jq)] (q G N*), (28) 



F p , q := E 

when p goes to infinity. This step relies on Theorem 2.1 in [GLOlj . which is a limit theorem 
for random walks on an affine group and generalizes theorems A and B in |LPP97j . 
Finally, we prove that the limit of F^ qt ^ q , when q — > 00, and a^(t) both have the same 
asymptotic behavior when t goes to infinity. 

Remark 1 . Let us mention limit behaviors in the discrete setting, which have comparable 
forms as our results. A BPRE (X n ,n £ N) is an integer valued branching process, speci- 
fied by a sequence of generating functions (/ n , n G N). Conditionally on the environment, 
individuals reproduce independently of each other and the offsprings of an individual at 
generation n has generating function f n . We present the results of Theorem 1.1 in jGKOO] 
and Theorems 1.1, 1.2 and 1.3 in |GKV03] (to lighten the presentation, we do not specify 
here the moment conditions). 

In the subcritical case, i.e. E[log(/Q(l))] < 0, we have the following three asymptotic 
regimes when n goes to 00, 

P(X n > 0) ~ ca n , as n — > 00, 
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for some positive constant c and 



a n = E 



/o(l) 



a n = n-^E f' {\) and a n = n" 3 / 2 min E (^(l)) 6 
L J \0<s<l L 



if respectively E[/g(l) log(/o(l))] is negative, zero or positive. 
In the critical case, i.e. E[log(/Q(l))] = 0, we have 



P(X n > 0) ~ cn" 1/2 , 
for some positive constant c. 



as n — > oo, 



5.1 Discretization of the Levy process 

The following result, which is a direct consequence from the definition of Levy pro- 
cesses, allows us to concentrate our attention on A Ptq , which was defined in ([2] 



Lemma 9. Let t > 1 and g 6 N*. There exist two random variables C^ qt ^ q and D^ qt j_i q 
such that 

ie^H/^ LgtJ _ 1)g <^ e -^ ds <ie^M/« q?tJ)? , 



and for every (p, q) £ N x N* , 



D 



(d) 



p,q 



U\/ q A p ,q and C Pj(] — V 1 ^A Pj(j , 



where the couple of r.v. (Ui/ q ,Vi/ q ) is independent of A Ptq and has the same distribution 
as (e 1 /i , e va). 



Proof. Let (p, q) be in N x N* and s in 



K s < K p/q + \i\/q+[a 
Now introduce 



,(+)l 



{p+1)/q ° p/q \ and K s >K p/q -Wq-[al p l 1)/q -a^]. (29) 



A" 



(i) 



Viz - K p/q + l a (p+l)/q u P /qi u l, 

Then, we have for all (p, q) G N x N 



(+) 



a^h-aiV and = K n 



(-) 



p/9 - iY P/9 Hp+1)/<? °p/«/ J ^ 1/9 ' 



(2) 1 



Moreover, tf« ) is independent of and (A< 2 \ A^ , A^) is hide- 

pendent of Finally, the inequalities in (|29|) lead to, 



L«*J— i 



l e ^(| 7 | /g+<7 (+)) e -^W < /■ < 1 ( 



i=0 



which ends the proof. 



□ 



18 



5.2 Asymptotical behavior of the discretized process 



We first recall Theorem 2.1 in |GL01j . in the case where test functions do not vanish. 
This result is the key to obtain the asymptotical behavior of the discretized process. 

Theorem 10 (Giuvarc'h, Liu 01). Let (a n ,b n ) n >o be a (M?j_) 2 valued sequence of iid 
random variables such that E[log(ao)] = 0. Assume that bo/(l — ao) is not a.s. constant 
and define Aq = 1, A n = YYkZo a k an d B n = ^2^Zo -^i^ii for n > 1. Let n, k, A be three 
positive numbers such that k < A, and 4> and ip be two positive continuous functions on 
R+ such that they do not vanish and for a constant C > and for every a > 0, b > 0, 
b' > 0, we have 



(a) < Ca K , $(b) < 



C 



and \4>(b) -${b')\ < C\b-b'\ T >. 



(1 + 6) A ' 
Moreover, assume that 

E[a£]<oo, E[ao' ? ]<oo, E[b v ] < oo and E[a~ v b~ v ] < oo. 
Then there exist two positive constants c((p, ip) and c(ip) such that 



lim n 3/2 E 



j>(A n )i>{B n ) =c(4>,i;) 



and 



lim n 1/2 E 



i>{B n ) = c$). 



Recall the definition of A Ptq and F P)Q in ()26p to ()28|) . The three following Lemmas 
study the asymptotic behavior of their expectations in the regimes (ii), (hi) and b/. 

Lemma 11. Assume that |<^-(0+)| < oo, 6 max > 1 and 4>'k(1) = 0. Then there exists a 
positive and finite constant C2(q) such that, 



F p ^C F c 2 (q)( P /qr 1/2 e^'<W, as 



p^oo, 



(30) 



and 



E 



(A p>q /q)-W\ ~ 02 (g)(p/g)- 1 /2 e (p/9)^(i) j as p ^oo. (31) 
Proof. Let introduce the exponential change of measure known as the Escheer transform 
dP (A) 

Ft 



dP 



for a e [o,e max [, 



(32) 



where (J-t)t>o is the natural filtration generated by K which is naturally completed. 
The following equality in law 



Ap^q — e 



V V 

P(K p/q -K i/q )\ (Jp-^fVp^ 



i=0 



i=0 



leads to e -(p/9)^(i) E Ap^P] = E« \A~} q /P , where A p>q = Yl=^ Kl/q - Let e > such 
that holds and observe that A Pjq > 1 for every (p, (/) in N x N*. Thus, 



E (i) < E (i) A-i/P < E (i) 



inf e 

ie[o,p]nN 
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Since 0^(1) = and E[ifj / ] < oo, Theorem A in |Koz76] yields 



E (i) 



inf e Kl / q 
ie[o,p]nN 



C q (p/q)- 1/2 , 



as p — > oo, 



where C„ is a positive finite constant. Define for z > 1 



Note that there exists po £ N such that for p > po, D q (l,p) < 2C q . 

Our aim is to prove that D q (l,p) converges to a finite positive constant ^2(9)- Then, 
we introduce an arbitrary x S (0, (Cf /M) 1//£ q~ 1 /^) and apply Theorem If 01 with 

$(z)=F(z), 4>(z) = z 1 /W\ k, A) = (f , 1/ (2/3), f / /3). 

Under the probability measure P^, (a n ,^n)n>o = (e /3( - K ^ n+i y''~ Kn ^\ x~Pq~ 1 ) n >a, is an 
iid sequence of random variables such that E^ 1 ) [log(ao)] = 0, as (j)' K (l) = 0. A simple 
computation gives us 

This implies that the moment conditions of Theorem [TU] are satisfied. In this case, 



E^Wl = e (^(i-«-^(i))/ g < ^ 



n-l 



B n = q- 1 x-PJ2 e$Ki/ "> neN*. 
i=0 

Thus, there exists a positive real number b(q, x) such that 

(p/<z) 1/2 E« |>(x-%,,/g)l -> 6(g, x), as p -> 00. 



Now, we define limin^^oo D q (l, n) = D„ and limsup n _ i>00 D q (l, n) = D q . Taking expec- 
tation in (|37|) yields 

{p/q) 1 / 2 ^ ^(x-PAfJq^ -C F xq l ^D q {l,p)\ < Mx 1+£ q( 1+£ ^ D q (l + e,p). 

Let n& and be two increasing subsequences in N such that, 

D(l,n k ) — > D q and D q (l,m k ) — > D q . 

k—^oo k—>oo 

As D q (z,p) is decreasing with respect to z, we have for all k in N, 

(C F xqW + Mx 1+e q<> 1+e W)D q (l t n k ) > (n k /q) 1 ' 2 ^ [F(x^A nk Jt 

and 

(CpxqW - Mx 1+£ q( 1+e »P)D q (l,m k ) < {m k / q ) 1 ' 2 ^ A mk Jq) 
This implies, taking k — > 00, that 



> 0, A, < 
20 



9 - Cfxq 1 //? - Mx 1 + £ g( 1 + £ )// ? 



< 00, 



and 



D q -D q < 



4MC (? xV 1+£)//3 
C~f ' 



Finally, letting x — > 0, we get that D q (l,p) converges to a finite positive constant d2(q), 
which gives (f3Tj) . 

Using ([37]), we get 



E 



F m - C F (A P}q /q) 



-1//3 



< E 



(A p Jq)-^ W 



so (I30p will be proved as soon as 



E 



PA 



o E 



as p — > oo. 



From the Escheer transform (|32p . with A = 1 + e, and the independent increments of 
K, we have 



E 



A-(l+e)/f} 

p,q 



e {p/q)4>n (i) E (i) [(Ve - ^*/* 



-e/0 



< e (p/<l)<t>K 
= e (p/9)te(l) E (l) 



-' (1) E (1) f inf e e ^/9 inf 

Lo<i<Lp/3j L2p/3J<j<p 



^ e /3(^p/ 9 -^/g) 

{K p/q -K j/q ) 



inf e ei ^/9 

0<«< Lp/3J 



E (l) 



inf e ^'/a 

0<j<[p/3j 



Using (|25p . we observe that E^(K 1/g ) = and E^(Kf /q ) < oo. We can then apply 
Theorem A in |Koz76j to the random walks (— iQ/ g ).;>i and (eKi/q)^. Therefore, there 
exists c > such that 



E 



^ 1+e)//3 ] < (c/p)eWto*M = o (E [a~^]) , 



as p — > oo. 



Taking C2(q) = d,2{q)q 1 '^ leads to the result. 



□ 



Remark 2. In the particular case when f3 = 1, it is enough to apply Theorem 1.2 in 
|GKV03| to a geometric BPRE (X n ,n > 0) whose p.g.f's satisfy 



fn(s) = J^Pn^f 



Pn 



k=0 



1 - q n s ' 



with l/p n = 1 + exp {/3 {K {n+1)/q - K n/q )}, q n = l- p n . Using E(ArJ) = P(X P > 0) 
and log/^l) = Ki/ q , allows to get the asymptotic behavior of E,(A~V) from the speed 
of extinction of BPRE in the case of geometric reproduction law (but we need the extra 
assumption 4>k(2) < oo). 

Recall that r is the root of 4>' K on ]0, 1[, i.e. 4>k{t) = mino< s <i<^.K-(s). 

Lemma 12. Assume that 0^(0+) < 0, 4>' K {1) > and /3 + 1 < max . Then there exist 
two positive and finite constants d(q) and c%{q) such that 

~ c 3 (q)(p/qy 3/2 e^^ KiT \ as p -> oo 



P.9 



and 



E 



(A^/g)" 1 ^] ~ d{q)(p/q)~ z / 2 e^ lq)MT \ as p -»• oo. 



(33) 
(34) 
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Proof. To prove these asymptotic behaviors, we apply Theorem 1101 with. 

$(z)=F(z), 0(z) = z T /P z>0, (n,K,\) = (l,T/P,l/P)- 

Under the probability measure F^ T \ {a n ,b n ) n >o = (e~ l3 ^ K< - n+1 ^ q ~ Kn ^),q~ 1 ) n >o, is an 
iid sequence of random variables such that E( r )[log(ao)] = 0, as 4>' k (t) = 0. The moment 
conditions 



e -*K{r)/q < and E (r) = e &K(f>+r)-*K(T))/q 



enable us to apply Theorem 1101 In this case, 

n-1 



< OO, 



i=0 



Then there exists 03(5) > such that 

E [F{ApJq)] e ~(p/9)^W = eW [F(^ Pi9 /g) e - r ^/«] ~ c 3 (g)(p/g)- 3 / 2 , as p -»• 00. 



This gives (|33j) . 

In order to prove 



E 



as p — y 00 



for d(q) > 0, we follow the same arguments as those used in the proof of Lemma [TT1 In 
other words, we define for z > 1, 

D q (z,p) = (p/qf^e-W^x^E \a~^ 

which is decreasing with respect to z. We obtain the same type of inequalities as in the 
Lemma [11] for the random variable A instead of A. 

Again we take e > such that (j47j) holds. Then from Lemma 7 in [Hir98], we know 
that there exists C q > such that for p large enough, 



E 



Pi 9 



< E 



A- 1 '? 



p,q 



< E 



inf e Ki l q 

«G[0,p]nN 



then we use Theorem 1101 to get <C lim inf n— >oo 

D q (l,n) = limsup n _ +00 L'g(l,n) < 00, 
which ends the proof. □ 

Lemma 13. Assume that (p' K (0+) = and (3 < 9 max . Then there exist two positive and 
finite constants b(q) and 04(g) such that 



and 



E 



F Piq ~ c A (q)(p/q) 1/2 , as p ->■ 00, 



(Apjq)- 1 ^} ~ 6(g)(p/g)- 1/2 , as p ^ 00. 



(35) 
(36) 
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Proof. The proof is almost the same as for the Lemma [T2l We first apply Theorem 
[TU1 to the same function -0 and sequence (a n ,6 n ,) n >o defined in Lemma [T2l but with the 
probability P instead of p( T ) . Then, we get 



E 



as p — > oo. 



We then define for z > 1, 



D q (z,p) = (p/q) 1 / 2 E\A-^ 



and from Theorem A in |Koz76j and Theorem I1UI . we obtain that D q (l,p) has a positive 
finite limit when p goes to infinity. □ 

5.3 From the discretized process to the continuous process 

Up to now, the asymptotic behaviors of the processes depend on the step size 1/q. 
By letting q tend to infinity, we obtain our results in continuous time. Recalling the 
notations (|30p to (|36p . we prove the following limits : 

Lemma 14. There exist five positive finite constants b, d, 02, C3 and C4 such that 

(b(q),d{q),c 2 (q),c 3 (q),C4(q)) — > (b, d, c 2 , c 3 , c 4 ), as q ->■ 00. (37) 
Proof. Let us first prove the convergence of d(q). From Lemma El we know that for every 



-(i)-|tI)/9]e 



A n q,q/q 



E 



n 



-3/2f,n<j> K (T) 



< 



< 



j^e-^du 



-VP- 



n -3/2 e n<f> K (T) 



,(<P+(l) + h\)/q E 



A 



nq—l,q 



In 



-VP 



n -3/2 e n</) K (T) 

A direct application of the Lemma [22] with 

a(q)=d(q), c~(q) = e^-W-W)/^ an d c+fa) = e^+W+WV?, 

yields that d(q) converges as q — > 00. 

Similar arguments lead to the convergence of b(q) and we now prove the convergence 
of C2(q), cs(q) and 04(g). Again the proofs of the three cases are very similar, so we only 
prove the second one. From Lemmas [9l and [T2| we know that for every (n, q) £ (N*) 2 , 



E 



F (e^vf /q A nq Jq)] < a F (n) < E ^(e^^U^A^Jq 



V1 1 



Using (|50p and dividing by n 3 ^ 2 exp(n^ft-(r)), we obtain 



F nM + ME 



-h\/qy 



V<i 



E 



-3/2 p n<t> K (r) 



< 



aF\n) 



-3/2 p n<j> K (T) 



n 



F, 



nq—l,q 



+ ME 



< 



Vi 



E 



An 



9-1.9 



n -3/2 e n<f) K (T) 
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Lemmas [T2l 1221 and equation ([50|) . where a(q) = 03(9), 

c"(g) = 1 ^ — '- and c + (<?) = H - 



C3(g) 



C3(g) 



yield the result. 



□ 



5.4 Proof of Proposition [8] 



Proof of Proposition a/ (i). Recall from Lemma II. 2 in [BLGOOj that the process {Kt ■ 
Ku_ s \-,0 < s <t) has the same law as (K s ,0 < s <t). Then 



e -P K °ds= I 
Jo 



*ds. 



We first note that for every q S N* and t > 2/q, Lemma [9] leads to 



E 



t \ -VP 



< E 



2/q 



-1/0 



-? K °ds 



< q^e^E^^A^ 

1//3 f <Ml) + |7| +<&(!) -, 
= q ,H exp I ) < 00, 



where <pfc was defined in (f23|) . Hence using the change of measure ([32]) . with A = 1, we 
have 



E 



t \ -W 



E 



t \ -VP' 
e? K -d8 







3 i0K(l) E (l) 



t \ -w 

ef* K 'ds 



The above identity implies that the decreasing function 1 1— > E^[(/q e /3Ks ds)~ 1 ^ 13 ] is finite 
for all t > 0. So it converges to a non negative and finite limit ci, as t increases. This 
limit is positive, since under the probability P^ 1 ) , K is still a Levy process with negative 
mean E^-\K\) = 4>' K (1) and according to Theorem 1 in |BY05j . we have 



e 0Ka ds < 00, 



(i) 



-a.s. 



It remains to prove that 



a F (t) ~ C F E 



ds 



-1//3 



as t — > 00. 



Recall that 9 max > 1 and <j>' K (l) < 0. So we can chose e > such that (|47|) holds, 
1 + e < 9 max: (j) K (l + e) < 0jf (1) and <f/ K (l + e) < 0. Then 



/ 6 



Thus, we just need to show that 
E 



* \ -1/j9 







t N -(l+e)/^ 



o(e*^ (1) ), as i^oo. 



24 



It is achieved by a new change of measure (|32p . with A = 1 + e, 



E 



t \ -(l+e)/j8- 



E 



pi \ -Q-+e)/P 

( I e? K >d8 



3 t0 K (l+e) ]E (l+£) 



o 



Again using Lemma [9j we obtain for t > q/2, 
J e~ p ^ds 



E 



-/3K S 



, (1+6)7/9 ^x(l + g) + |7l(l + £) + ^(l + e) \ 



which ensures that the decreasing function i i— )• E( 1+£ )[(J"q e pKa ds) ( 1+e )/^] is finite for all 
i > and gives the result. □ 

Remark 3. In the particular case when j3 = 1, it is enough to apply Theorem 1.1 in 
[GKV03| to the geometric BPRE (X n ,n > 0) defined in Remark [2] to obtain the result. 

Proof of Proposition^ a/ (ii), (Hi), and b/. Proofs are similar for the different regimes, 
so we only focus in the proof of the regime in a/ (iii). 

Let e > and q £ N* such that q > 1/e and (1 - e)c 3 < c 3 (q) < (1 + e)c 3 . Then for 
every t > 1, 



< F L9tJ _ 1|9 +E[F(D LgtJ _ 1)ff e-^W/V?) 

Applying ([50jl . we obtain : 



LgtJ-l.g- 



E 



'[.9*1,9 



\qt-l\,q 



< (i _ e -(l7l-^(i)))ME[( J 4 L(?tJig /g)- 1 //3 

< (e e ^l+*+W) - l)ME[(A LgtJ _ l!g / g )- 1 ^ 



When i goes to infinity, we can bound both terms by 



h ( £ ) t -W e t4>K(T) 



2Md(e £ ^ +< ^ + ^ - e~ £ ^~' t> -^)e~ E ^ K(:r ^ 



3/2J<t> K (r) 



(38) 



where and are defined in (|23p . and /i(e) goes to with e. On the other hand, for 
t large enough 



(1 - 2e)c 3 t- 3 / 2 e t ^ < F lqtu < F [qt] _ hq < (1 + 2s)c 3 t~ 3 / 2 e^^, 
which ends the proof of Proposition [8l 



□ 



6 Auxiliary results 

This section is devoted to the technical results which are necessary for the previous 
proofs. 
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6.1 A local martingale 

To prove the expression of the Laplace exponent of Z in Theorem [TJ we used a local 
martingale, which is determined below via Ito's formulae. 

Lemma 15. For every F £ C 1 ' 2 (R+,1R), conditionally on A, 

f't £l2 pt 

F(t,Z t )-F(P,Xo)- J F(s, Z s )a 2 e- 9S - A °Z s ds - J -F(s,Z s )ds 

- £ J™Y s (F(s,Z s + ze-f s " A - F(s,Z s ) - ^-F(s,Z s )ze-^- A ^fi(dz)ds, 
is a local martingale. 

Proof. First, we apply Ito's formula to Z% and 

Z t = Y + f e-° s - A W2o-*Y s dB s + g f Ye^^ds 
Jo Jo 

+ f (°° I 8 e- 9S - A °-zN (ds,dz,du) -g f Y s e~ 9S - As ds 
Jo Jo Jo Jo 

+ [ [ (ze- 9S - As — logW -e- 9S - A ^)y s _^i(d5,dz) 

JO J[0,oo) V ' 

= Y + e~ 9S - A ^2a 2 Y s dB s + / / / e" 9 ^- zN (ds, dz, du), 
Jo Jo Jo Jo 



which is a local martingale. 

Now, again from Ito's formula, we have 



F(t, Z t ) = F(0, Y ) + J* ^F(s, Z s )e- 9S - A * V^Y s dB s 

pt £\2i pt o 

+ J ^F(s,Z s )o- 2 e- 9S - A °Z s ds + J -F(s,Z s )ds 

pt f OO p "¥ 

+ J J J ' [F(s,Z s ^ + ze~ 9S - A °-) - F(s,Z s ^N (ds,dz,du) 

pt p OO e ~\ 

+ J J Y s (f(s,Z s + ze~ 9S - A °) - F(s,Z s ) - — F{s,Z s )ze~ 9S ~ A ^^(dz)ds. 
It yields the expected local martingale, conditionnaly on A, and ends up the proof. □ 



6.2 Existence and uniqueness of the backward ODE 

The Laplace exponent of Z in Theorem [1] is the solution of a backward ODE. The 
existence and uniqueness of this latter are stated and proved below. 

Proposition 16. Let 5 be a cadlag function on M + of bounded variation. Then the 
backward ODE admits a unique solution. 

The proof relies on a classical approximation of the solution of ((HJ) and the Cauchy- 
Lipschitz theorem. When there is no accumulation of jumps, the latter ensures the 
existence and uniqueness of the solution between two successive jump times of 5. The 
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problem remains on the times where accumulation of jumps occurs. Let us define the 
family of functions 5 n by deleting the small jumps of 5, 



5 t = St - £ (s. 



s<t 



s- 1 *-{\5 s -5 s -\<l/n}- 



We note that i/jq is continuous, and s i— > e 9S+<5 " is piecewise C 1 (M + ) (with a finite 
number of discontinuities). From the Cauchy-Lipschitz Theorem, for every n G N* we 
can define a solution vf(.,X, 8) continuous with cadlag first derivative of the backward 
differential equation: 



d 

dS v * 



n {s,X,5) =e 9S+5 °ip (e- 9S - 5 "v?{s,\,8)), 0<s<t, v^(t,X,S) = X. 



We want to show that the sequence (vf(., X, <5)) n >i converges to a function vt(.,X,S) 
solution of KSJ. 

Then the proof of this proposition follows from the next result. Let t be fixed and 

S:= sup \e 9S+s ",e~ 9S ~ s °}. (39) 

s€[0,t],neN* J 

Lemma 17. For every A > 0, 
(i) we have 



inf v?(s,X,5)>0, (40) 

0<s<t, Kn 



(ii) there exists a positive finite constant C such that for all IS 1 < rj < k < XS, 

0<M*)-Mv)<C(K-ri). (41) 

Proof. First, we observe that S is finite. Now, using that x \— > e~ x + x is increasing on 
and the Taylor-Lagrange's formula to x t— > e~ x , we can check that for all < rj < k 
and x > 

< e ~ KX ~ e ~ VX + ( K ~ < l + K ~ v (42) 
(k — rf)x 2rje 

In order to prove (i), we first note that ipo(X) > for all A > 0. Moreover, if there 
exists < so < t, such that v^(sq, A, 5) = 0, then vf(s, X, S) equals for every s G [so,t] 
since V'o(O) = 0. Hence vf(., X,5) is non decreasing and 

<(s, A, 5) £ (0, A] for all s G [0, t] and n > 1. (43) 

Moreover -00 is increasing and for all s G [0,t], and re > 1: 



e 



where S is the constant defined in (|39|) . On the other hand, note that for all x > 0, 

< e~ x - 1 + x < x Ax 2 , (44) 

then for all < v < AS 1 , we get 



l 

re := Bv. 



Tpo{v) < XS^ J x 2 n(dx) + cr 2 ^ + y xn(dx 
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Putting all the pieces together, we deduce that for all < s < t, 



and since vf(t, A, 5) = A, we have 

v?(s,X,6) > Ae 5 - 8 ^-*), for all s G [0,t]. 
We get (i) by defining I = Xe~ SBt . 

Finally, we note that for all IS 1-1 < r] < k < XS, 
ipo(n)- tpo(v) 

= a 2 (K 2 - V 2 ) + (e~ KX - e - "* + (k - rj)xj fi(dx) 



<<r(Kr-rf) + {,K-rj)[\ + 



2r]e 



xfi(dx) 



< 



l 

+ (k — 77)77 / x 2 fi(dx) + (k — r/) 2 / x 2 fjL{dx) 
Jo Jo 

2A5tJ 2 + 2A5 / zVdc) + (1 + — ) / x/i((fx) (ft-?/), 
Jo 2/e 7i J 



which proves part (ii). 

We can now prove the result of existence and uniqueness. 



□ 



Proof of Proposition [TR We now prove that (vf(s, A, S), s £ [0, i]) n >o is a Cauchy se- 
quence. For sake of simplicity, in what follows we denote v n (s) = v^(s, X,5), and for all 
v > 0: 

i/> n (s,v) = e 9S+s "i;o{e~ gs ~ 5 "v) and 4>°°(s,v) = e 9S+<5s Vo(e~ 9S ~^) ■ 
We have for any < s < t and m,n > 1: 
\ v n ( s )-v m (s)\ 

= j iP n (u,v n (u))du- j i) m (u,v m (u))du 

J S J S 

< I R n {u)du+ ! R m {u)du+ I ijj co {u,v n {u))-^ 0O {u,v m {u)) du, (45) 

Js Js Js 

where for any u G [0,t], 

R n (u) = ^ n (u,v n (u))-^ co (u,v n (u)) 



(«)) 



e 5 " - e 5 ° 
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Moreover, from ([39]) . (HPT) , and (141 p. we obtain 



< S'CA 



+ e 9 Vo(A5) 



e 5 ? - e *« 





e -5? _ e Ss 




e 5 " - e s ° 


) sup { 




1 ae[0,t] L 









which implies 



sup//" R n {u)du,s G [0,t]\ < ts n -> 0. 
Using similar arguments as above, we get from (|4ip . 



Prom (|45p . we use GronwalPs Lemma (see Lemma 3.2 in |Dyn91| or Lemma 4.6 in |BS47] ) 
with 



Rm,n( s )= I R n {u)du+ f R m {u)du, 

J s J s 



to deduce that for all < s < i, 

\v n (s) - v m (s)\ < R m , n (s) + CS 2 e cs2 ^ 
Hence for every no G N*, 



sup \v n (s)-v m {s)\ <t 

m,n>no,s£[0,t] 



1 + CS 2 e csH t 



sup (s n + s m ). 

m,n>no 



Thus (v n (s),s G [0,t]) n >o is a Cauchy sequence under the uniform norm. Then there 
exists v a continuous function on [0, t] such that v n — >■ as n goes to oo. Now we prove 
that v is solution of the equation ([8]) . By continuity, v satisfies (|39p and (|41[) . Then for 
any s G [0, t] and n G N*: 



v(s) - / ^°°(s,t>(s))ds-A 



< 



v(s)-v n (s) + f ^{s,v(s))-^ n (s,v{s)) ds+ f i; n (s,v(s))-^ n (s,v n (s)) 

J s J s 



ds 



< ts n + (1 + CS 2 ) sup { v(s) -v n (s) , s G [0,t]}. 
so that letting n — > oo yields 

«(«)- / ip°°(s,v(s))ds — A 



0. 



It proves that f is solution of ([5]), whereas the uniqueness follows from Gronwall's lemma. 

□ 
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6.3 An upper bound for i/jq 

The study of the Laplace exponent of Z in Corollary [2] requires a fine control of the 
branching mechanism ^o- 

Lemma 18. Assume that the process (gs + A s ,s > 0) drifts to +oo a.s. There exists a 
non negative increasing function h on M + such that for every A > 



^o(A) < \h{\) and J /i(V (9 * +At) ) dt < oo. 
Proof. Inequality (|44p implies that for every A > 0, 

POO 

■0o(A) < o- 2 A 2 + / (A 2 x 2 l {Aa .< 1} + Xxl {xX>1} ) n(dx) 
Jo 

< (a 2 + J x 2 fi(dx)^j A 2 + A 2 1{a<i} x 2 fi(dx) + \J^ xfi(dx). 
Now, using condition (|13p we obtain the existence of a positive constant c such that 



A / x(i(dx) < A log(l + l/A)~ (1+e) / x log(l + x) l+t {i{dx) < cA log(l + l/A)~( 1+e ) 
Ji/\ Ji/\ 



Next, let us introduce the following function /: 

, log(l + x) 1+e , F1 . 

f(x) = — — , forxG l,oo). 

x 

By differentiation, we check that there exists a positive real number A such that / is 
decreasing on [A, oo). Therefore, for every A < 1/A, 

l/A ,1/A 

X 2 x 2 fi(dx) < A 2 / x 2 fi{dx) 
J A 



< Alog(l + l/A)-W/(l/A) f 'y ^) 

/•l/A 

< Alog(l + l/A)- (1+e) / xlog(l + x) 1+£ //(cix) 

J A 



Adding that for A G [1,1/A], fi /X \ 2 x 2 n(dx) < A 2 J^x 2 fi(dx) < \ 2 Af™xfi(dx) and 
using again condition ()13|) . we deduce that there exists a positive constant d such that 
for every A > 0, 

V'o(A)<c / (A 2 + Alog(l + l/A)-( 1+e )). 

Since A 2 is negligible with respect to Alog(l + l/A)~( 1+£ ) when A is near or infinity, we 
conclude that there exists a positive constant c" such that 

^o(A) <d"Alog(l + l/A)-( 1+e >. 
Defining the function h(x) = c"xlog(l + 1/ x)~^ 1+£ > , for x > 0, we get: 

f°° / \ f°° 

0< J h^ 9t+At ^dt<c" J (gt + A t y^dt, 
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which is finite since the process (gs + A s , s > 0) goes linearly to +00 a.s. More precisely, 
one can find e > such that (gs+A s — es : s > 0) has positive expectation for s = 1, which 
ensures that it goes to 00 a.s. and there exists L > —00 a.s. such that gs + As > L + es 
a.s. This ends the proof. □ 

6.4 Extinction versus explosion 

Let us here check that Yj can be properly renormalized as t —> 00 on the non-extinction 
event. 

Lemma 19. Let Y be a non negative Markov process which satisfies the branching prop- 
erty. 

Assume also that there exists a positive function a t such that for every xq > 0, there 
exists a non negative finite random variable W such that 

a t Y t >W a.s, P X0 (W > 0) > 0, a t 0. 

t— >oo 

Then 

{W = 0} = {Y t >0\ ¥ X0 a.s. 

Proof. First, we prove that 

P xo (limsup Y t = 00 I lim sup Y t > 0) = 1. (46) 

t— >oo t— >oo 

Let < x < x < A be fixed. As a t ->■ and ¥ X (W > 0) > 0, there exists t > 
such that a := ¥ x (Yt > A) > 0. By the branching property, the process is stochastically 
monotone as a function of its initial value. Thus, for every y > x (including y = xo), 

Wy(Y to >A)>a>0. 

Let us define the following stopping times 

T := 0, T i+1 = mf{t > T t + t Q : Y t > x} (i > 0) 

For any i 6 N*, by strong Markov property 

^x (Y Ti +to >A\{Y t :t< Ti), Ti < 00) > a. 

Conditionally on {lim sup t ^. 0O Yi > x}, the stopping times Tj are finite a.s. and for all 

< x < x < A, 

IPx (Vi > : Y Ti+t0 < A, lim sup Y t > x) = 0. 

t— >oo 

Then, ¥ xo (lim sup^^ Yt < 00, lim sup^^ Yt > x) = 0. Adding that {lim sup^^ Y t > 
0} = U a . e(0i2 , ]{limsup t ^ oo Y 4 > x} yields ([MD- 

Let us now consider the stopping times T n = inf{£ > : Yt > n}. We get by strong 
Markov property and branching property, 

¥ X0 (W = 0;T n < 00) = E X0 (l Tn<oo F Y (W = 0)) < F n (a t Y t — > 0) = Pi(a 4 Y t — )■ 0) n , 

which goes to zero as n — > 00, since Pi {a^Yt 0) = Pi(W = 0) < 1. Then, 
= F X0 (W = 0;Vn : T n < 00) = F X0 (W = 0,limsupY" t = 00) = ¥ X0 (W = 0,limsupY t > 0), 

t— >oo t— >oo 

where the last identity comes from (|46p . It completes the proof. □ 
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6.5 A Central limit theorem 

Finally, we need the following central limit theorem of Levy processes in Corollary [3J 
Lemma 20. Under the assumption (f 14\ ) we have 

py/t t^cv 

Proof. For simplicity, let r\ be the image measure of v under the mapping x h- > e x . The 
assumption (fbij) is equivalent to Jj a .| >1 x 2 n(dx) < oo, or EfA 2 ] < oo. 

Next, we define T(x) = ??((— oo,— x)) +r/((x,oo)) and f7(x) = 2 J* yT(y)dy, and 
assume that T(x) > for all x > 0. According to Theorem 3.5 in Doney and Mailer 
[DM02] there exist two functions a(t), b(t) > such that 

gt + A t -a(t) d Aj(n U(x) 

—- > N{0, 1), if and only if > oo. 

b(t) t^oo x z l (x) x->oo 

If the above condition is satisfied, then b is regularly varying with index 1/2 and it may 
be chosen to be strictly increasing to oo as t — > oo. Moreover b 2 (t) = tU(b(t)) and 
a{t) = tA(b(t)), where 

A{x) = g+ f ZT?(dz)+?/((l,oo))-?/((-oo,-l))+ / (v{(y, oo)) -t/((-oo, -y)))dy. 
J{\z\<i} Ji v 7 

Note that under our assumption x 2 T(x) — > 0, as x — > oo. Moreover, note 

z 2 r){dx), 



and 



U(x) = x 2 T(x) + / z 2 r](dx) + / 

J(-x,0) J(0,x) 

[ » K d I ) + ,(,((,,=.))-,((-=^-x))). 



C/(x) 



'{|z|<a:} 

Hence assumption (|14p implies that 

J (-00,00) 

Therefore, we deduce [/(x)/(x 2 T(x)) — > oo as x — >• oo, ~ p\/i and a(t) ~ mi, as 
i — > oo. 

Now assume that T(x) = 0, for x large enough. Define 

'gt + At - a(ty 



*(A,t) = -logE 



exp < iA 



6(t) 

where the functions a(i) and 6(t) are defined as above. Hence, we can write 

*(\,t)=t[ (l-e^) x + ^-x + ^-x 2 )i 1 {dx) + t [ (l - eW) x ) v (dx) 

J{\x\<b(t)} v b{t) 2b CO y Aw^K*)} V 7 



x 2 r/(dx) + iAt ( rj(b(t), oo) — oo, —b(t)) 

{\x\<b(t)} 



2b 2 (t) 

Since T(x) = for all x large, 6(t) — > oo and t~ l b 2 {t) — >• p, as f — > oo, therefore 

t— »oo z 

which implies the result thanks to Levy's Theorem. □ 
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6.6 Two technical Lemma 

We now give two technical lemmas, useful in the proofs of Section 

Lemma 21. Assume that F is defined as in 1126}) . then there exist two positive constants 
i] and M such that for all (x, y) in and e in [0, rj\, 



F(x) - C F x 



< F{x) < M(x + l)~ 1//3 , 
F(x)-F(y) < M\x-y\, 



F(x) - F(y) 



< M 



X -1/P _ y-l/P 



(47) 
(48) 
(49) 

(50) 



Proof. We only prove (|50p . since the others inequalities follows from straightforward 
computations. We define the function h : x £ M + i— > (1 + x) 1_? /i(x) and let < x < y. 
Then, 



0<(F(x)-F(y))/C F < (x + 1) 



(y + l)-^)+(l +y ) 



-1//3-1 



+ 



/ i (x)|((l + x)- 1 / /3 - 1 -(l + y)- 



l/p-i 



h(x) - h(y) 
(51) 



1 + x 
x) 



l + y 

l + x 



i _ 2 l + y 



+ 1 



1 + y l + x 



As /3 £ (0, 1], we have the following inequalities : 

(i + x)- 1 //3-i_(i + y )-i//3-i < i i + y) -^-U( l + y\ l/ p 

< (i + y) 

< 2(x- 1 ^-y 

Moreover, the Mean Value Theorem applied to the function z £ M + i— >• (z + l)" 1 /' 3 on 
[x, y] ensures that 

X -{y + l)- 1 ^"!^ - x ) < ( x + I)" 1 // 3 - (y + I)" 1 /*. 
Now, we denote by the Lipschitz constant of h. The equation (|5ip finally gives 



0<F(x)-F(y) < C F (l + 2||/ i || 00 + A;/3)(x- 1 ^-^ 1 //3 



which ends up the proof. 



□ 



Lemma 22. Assume that the positive sequences (a n ,q)( n ,q)e^ 2 , fan q ) ( n ,q)eW an <^ (^n)neN 
satisfy for every (n,q) £ N 2 : 

anc? i/iai i/tere exisi too sequences (c~(g)) ge N ( c+ (</)<?eN swc/i t/iat 

lim a n g = c~(q)a(q), lim = c + (q)a(q), and lim c~(g) = lim c + (q) = 1 



q— >oo 



ITien i/iere exists a positive finite constant a such that 



lim a(gj = lim a'(q) = lim b n = a 

q— ¥00 q— >oo 71— >oo 
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Proof. Letting n go to infinity, we have for every q £N 

limsupfr n < c + (q)a(q) and c~(q)a(q) < liminf b n . 
Then letting q go to infinity, we obtain 

limsup6 n , < liminf a(q) and lim sup a ((/) < liminf b n , 
which ends the proof. □ 
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